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Abstract 

A systematic extension of the Monte Carlo (MC) algorithm, that solves the 
DGLAP equation, into the so-called the one-loop CCFM evolution is presented. 
Modifications are related to a ^-dependent coupling constant; transverse momentum 
dependence is added to the x-dependence of the parton distributions. The presented 
Markovian algorithm for one-loop CCFM evolution is the first step in extending it to 
other more sophisticated schemes beyond DGLAP. In particular, implementing the 
complete CCFM will be the next step. The presently implemented one-loop CCFM 
option will be a useful tool in testing the forthcoming MC solutions. Numerical 
results of the new MC are confronted with other non-MC numerical solutions. The 
agreement within the MC statistical error of ~ 0.1% is found. Also, numerical 
results for & T -dependent structure functions are presented. 
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1 Introduction 



The evolution equations (EVEQs) play a crucial role in the construction of any parton 
shower Monte Carlo program simulating the production of multiple gluons and quarks, 
in an approximate way, in perturbative QCD. In particular, EVEQs are in this context 
used to model the actual development of the tree-like cascade of primary emitted partons 
from the initial hadron. Such a cascade is then supplemented with an appropriate hard 
scattering matrix element and with the hadronization mechanism to form the complete 
parton shower MC code. 

The basic, and by far the best analyzed, EVEQs in QCD are the DGLAP equations [1]; 
originally formulated in the leading order, later extended to the next-to-leading order [2-5] 
and recently even to the next-to-next-to-leading order [6,7]. These equations perform 
systematic resummation in terms of the variable logQ, where Q is a hard scale, for not 
too small or too large values of the longitudinal momentum fraction x. 

Another important EVEQ is the BFKL equation [8]. Originally formulated in the 
leading order, it has been recently extended to the next-to-leading order [9]. Contrary 
to the DGLAP equations, it resums large log a; terms, therefore, describing the small x 
limit better than the DGLAP equations. In recent years attempts have been made to 
improve the DGLAP evolution equations in the small-x region by incorporating some of 
the BFKL features [10, 11]. However, these attempts have lead to rather complicated 
evolution equations. 

The third class of EVEQs in QCD are the CCFM equations [12]. They are formulated 
for the so called "unintegrated" parton distributions, which depend on the transverse 
momenta of partons k T in addition to x and Q. The main idea of the CCFM approach is 
to correctly describe not only the large-x region, where the summation of logg dominates, 
but also the region of small x, where the large logarithms ln(l/a;) are important. Thus, 
the CCFM equations effectively interpolate between the DGLAP and BFKL equations. 
One of the variants of the CCFM equations, called "one-loop" [13], is of special interest 
because of its simplicity and ease in extending to non-leading orders and to non-gluonic 
evolutions. 

The one-loop CCFM equation combines coherence effects (angular ordering of gluonic 
emissions) at large x with the transverse momentum ordering at small x [14]. In contrast, 
the "all-loop" CCFM formulation extends the angular ordering to the small x region. In 
its original formulation the CCFM scheme has been formulated only for gluonic cascades. 
In the "one-loop" approximation, as defined in ref. [15], it was extended to include also 
quarks into it. 

Finally, let us mention yet another class of equations, the IREE ones [16], which 
account for both the logarithms of Q and x by means of constructing and solving the 
so-called two-dimensional infrared evolution equations. 

In a series of recent papers [17-25], a new family of Monte Carlo algorithms that 
solve the EVEQs in QCD has been developed. The novelty of these algorithms is in their 
ability to include an energy constraint on top of the normal Markovian random walk-type 
evolution. This is of the utmost importance in the case of the MC simulations of narrow 
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resonance production processes initiated by initial state parton cascades. In parallel 
to the "constrained" algorithms, we have also developed the standard "unconstrained" 
ones, which are useful for tests and precise studies (for example, for fitting deep-inelastic 
structure function F2). In the very beginning of the above studies, the DGLAP evolution 
and some of its variants were the main object of the interest. However, from the previous 
discussion it is clear that it would be useful and interesting to extend them also to other 
types EVEQs, like those mentioned, having in mind the construction of the parton shower 
Monte Carlo at a later stage. 

In the presented paper we take a first step in that direction and show how the Marko- 
vian "unconstrained" DGLAP algorithm can be extended to the one-loop CCFM case. 
This extension involves two steps: the modification of the coupling constant which be- 
comes z-dependent and the inclusion of transverse momenta into the evolution. 

In the subsequent papers further extensions of the Markovian algorithms will be pre- 
sented, in particular to the all-loop CCFM scheme. All technical aspect of the imple- 
mentation of these "constrained" MC algorithms are discussed at length in ref. [26]. 
Summarizing, this work should be be treated as a warm-up exercise for the forthcoming 
more complete and also more sophisticated analysis. 

The paper is organized as follows. In Section 2 we recall some basic formulae on the 
DGLAP equations and describe in detail problems to be solved. Section 3 is devoted to 
the changes in the algorithm related to the modification of the argument of the coupling 
constant. Section 4 describes the extension of the algorithm to the /c T -dependent evolu- 
tion. In Section 5 some numerical results are presented, and finally, Section 6 contains 
summary and outlook. 

2 Framework 

In this paper we will heavily rely on the notation and formulas from [23] where the reader 
will find more details on the notation and description of the framework which we use. We 
start from the DGLAP equations^ 

1 

J x 

where x the energy fraction of the hadron carried by the parton of type K = gluon, quark, 
antiquark. The so-called evolution time is t = InQ = ln/i, with Q being large energy 
scale defined by the hard process probing parton distribution function (PDF) Dx{t,x). 

1 This is eq. (104) in ref. [23]. We shall provide explicit link to key formulas in this work in the 
following. 
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The iterative solution solution of the above equations rea> 



xD K (t, x) = e-* K ^xD K (t , x) + J2 dx o II / dU 9 ^ ~ 

n=l Q K ,...,K n ^i i=l 

n r -i n 



i=i 



8=1 



(2) 



where ii' = K n . In order to turn the series of eq. into the one-loop CCFM evolution 
we add two elements. First, the argument of coupling constant is made z-dependent: 
ct s (t) — ► Qj s (t+log(l— 2;)) and second, the evolution has to include the transverse momenta 
of emitted partons. As for the change of argument of a s there are several technical points 
to be clarified. 

(i) We restrict ourselves to the LL approximation only, so the kernels 7k 1 k^ 1 (U, zi) 
have the forrrjf] 

7 KJ (t,z) = 2P KJ ( Zl a s ), (3) 



Pkj(z, a s 



a { I>\t + \n(l-z)) 
2?r 



4%) 



(4) 



where functions P^j(z) are defined in eqs. [(A. 5)]. The factor of 2 in eq. ([3]) is related 
to the definition of the evolution time t = lng. The coupling constant a s (t) is defined as 
follows 

(ii) In the original papers [27-29] the shifted argument of a s has been used only for 
the diagonal kernels (Q— >Q or G— >G). We have, however, decided to apply the same shift 
for all kernels, including the Q<->G ones too. 

(iii) In the presence of the (1 — z) factor the argument of the coupling constant can 
become arbitrarily small. Therefore, to avoid the Landau singularity we impose the 
following infrared (IR) cut-off depending on the go parameter: 



9.(1 - Zi) > g - 



(6) 



It can be translated into 2j < 1 — g /gj in space of z- variables. The IR cut-off e = go/g; 
is not necessarily infinitesimal anymore, for close to go- 

(iv) As a consequence of a finite IR cut-off on the real emissions, the momentum sum 
rule 

j dzJ2^J K (t,z) = (7) 



2 See eq. (105) in ref. [23]. 

3 See eqs. (50) and (A.2) in ref. [23]. 
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gets violated (by ~ Q/qo corrections). In order to restore the sum rule we adjust appro- 
priately the virtual parts of the diagonal kernels 



? KJ (t, z) = -y 5 KK (t, g ) 6 KJ 6(1 -z) + ?%j(t, z), 
?%j{t,z)^? K j(t,z)Q{l-z-q e 



where 



1*kk= I dzzJ2?JK(t,z). (9) 

J T 



J 

In the following sections we will describe in detail all of the above technical points. 



3 The ^-dependent strong coupling 

In this section we will describe the MC implementation of the ^-dependent coupling 
constant. Let us recall the probability distributions of a single step forward in the Markov 
process. The Sudakov form factor is defined in the usual way 



M«-i5*o)= / dt'T s KK (t',t )= dt'J^ dzz?%(z) (10) 

and the probability of a step forward in z 

The probability of a step forward in t is given by the integral over dzi of the above duo 

du{t ^ K ^ = ^^.A^x^^-i^oJc-*^^"- 15 - . (12) 
In the above formulae we included also the new, finite cut-off Q). 

3.1 Singular part of the kernel 

Let us first calculate the Sudakov form factor only for the singular part of the kernels. 
The approximate kernels are now 

«; A,) = ^ " + ^ (1 " Z)) 9(1 ~ Z _~ z <t)] 4?« . K=G, tti , (13) 

where e(t) = exp(t — t). The virtual components will be discussed shortly. The kernels 
of eq. ([131) should appear in denominators of the expression for the compensating weight 



^Following eq. (69) of ref. [23]. 
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wp. In this particular case we are able to calculate the form factor analytically in an 

exact way 

ti 

*2 

= TT f dt' [ dlnx -4 ; 

/?o 7 J t'-lnA + lnx x ^ 

tl to-t' 

(14) 



|- / dt' pn(t' - In A ) - ln(t - In A )] = ^ I dt 1 (r' - r ) A 

po y po y 

T2 

f dr' e ^(r'-r ) = ^A 

J Po 

n 



(0) 



tl tl 

T'2 

^-A% [ dr' e *V - ro) = I (r' - ro - 1) e 

n 

2 



T2 



where 



A) 

g(r)= e '(r - to - 1) + e*» 

= (t - In A ){ln((t - In A )/(t - In A )) - 1} + (t - In A ), 1 } 

q{tq) = and we have used as previously t = t(q) = lng, r(t) = ha(t — lnA ), while 
t = r(to) = r(t(g )) is related to go in the IR cut 1 — z-i < e(qi, go) = qo/qi- The above is 
the part (dominant one) of the virtual form factor, but it suggests clearly another change 
of the evolution variable r — > g at the MC generation level of primary distributions, before 
applying the MC weight, 

Q(r(q)) = Q(q; 9b, A ) = e T <«> [r(g) - r(g ) - 1] + e T <»> , (16) 
which could increase the MC efficiency due to 

= dr e'(r - 7*) = df f\ z ^pM ^ ~ z ~ ^ (17) 

yo af(^) i-* 

The function £>(t) cannot be inverted analytically. However, we have implemented in 
the MC a simple and fast subprogram for the numerical inversion t(p) = £> _1 (t). Thus, 
the ^-function can be used to generate the variable t for every MC event. The efficiency 
of the MC is improved even more by performing importance sampling for the variable z 
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as well. To this end we calculate analytically the integral: 



ai 0) (t + ln(l- z)) 0(1- z - e (t)) 



dz 



71 



1 - Z 







— / d In x 



— ln(lnx + t — In Ac 



t — In A + In x Po 
where 

C(z) = ln(ln(l -z) + t-lnA ). 
The function ((z) can be inverted analytically, giving rise to 

z(C) = C -1 (z) = l-exp(e c -t + lnA ). 



lnz=0 



In x=to—t 



(19) 



(20) 



Therefore, the variable £ can be used to generate z according to the primary MC distri- 
bution. 

In order to check how this works in practice, we started from pure gluon-strahlung 
in the LL approximation. From the gluon-gluon momentum kernel we retained only the 
part of zy singular in (1 — z), which we denote as zy GG (z, t; t ), where 



a< 0) (* + ln(l-*)) A(0) 



A 



71 



GG 



1 1 
+ - 



Z Z 



9(1 -*-c(t)). 



(21) 



Note, that eq. (12T1) is given for D 5 , so also the 1/z singularity is present^]. Let us remark 
that the gluon momentum sum rule 



dz z'J ) GG (z 1 1; to) 







(22) 



is fulfilled by the form factor of eq. (THl) . As a result, the MC weight from this part of the 
algorithm is exactly equal to 1. In fact, the total event weight may slightly differ from 
1, if generation of x and K according the initial distribution D Ko (q ,x ) is done using 
FOAM [31] in the mode of weighted MC events. 



3.2 Non-singular terms in the gluon kernel 

In the previous subsection we have considered a simplified case of gluon-strahlung, retain- 
ing only the singular terms of the gluon kernel. Now we are going to extend this analysis 

'° This splitting function is just the kernel of the one-loop CCFM equation as formulated in Ref. [13]. 
The only difference is that in Ref. [13] the terms 1/(1 — z) and 1/z are multiplied by the running QCD 
coupling a s with different arguments, while in the above equation we use the same argument of a s for 
both terms. Such a choice is suggested by the NLL corrections within BFKL scheme, see e.g. [30]. 
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by including also the non-singular terms. It will be demonstrated on the example of the 
GG kernel. The case of the QQ kernel can be treated in the same way. The full LL 
gluon-gluon kernel corresponding to the real-gluon emission with ^-dependent a s reads 



7, 



a s (t,z) 



2C, 



71 



1 1 

1 — z z 



2 + z(l 



e{l-z-e{t)). (23) 



For the kernel corresponding to the gluon momentum distribution we get 



a s (t + ln(l - z)) 



2C f 



71 



+ F GG (z) 



e(l-z-e{t)), 



(24) 



F GG (z) = z[-2 + z(l - z)}. 



Enforcing the validity of the gluon momentum sum rule of Eq. (1221) . the following Sudakov 
form-factor is obtained 



$£(t 2) ti;t ) 



dt 



dz z7 GG (z, t; t ) 



(25) 



For the non-singular part of the gluon kernel, however, the integral over z cannot be 
calculated analytically. After integrating over t, for the non-singular part of the Sudakov 
form-factor exponent labeled as <&q, we obtain easily 



$ G , {t 2 ,t 1 ;t ) = -2C A 
Po 



l-exp(to-ii) 



dz In 



+ 



W 2Ca 



1— exp(to— *2) 

1— exp(to— ti) 
ti—to 



dz In 



due u In 



o 



t 2 - lnA + ln(l - z) 
ti - lnA + ln(l - z) 

t 2 - lnA + ln(l - z 

t - In A 
— In Aq — u " 



F GG (z) 
F GG (z) 



+ 



te—to 



due "In 



t\ — In A — u 
t 2 — In A — m 
t - lnA 



(26) 



GG\ 



1 GG{ 



where u = — ln(l — z). Numerical evaluation of the above one-dimensional integral can be 
done quite precisely (~ 0.01%) and rather quickly, as compared to the time of generating 
a single MC event, particularly in its latter form. The singular part, $q, of the Sudakov 
form-factor exponent is given by Eq. (1141) . The total Sudakov form- factor exponent is 
now 

$ G = + (27) 

The non-singular gluon kernel terms can be easily implemented in the forward Marko- 
vian algorithm of the previous subsection through appropriate MC weights. For the 
real-gluon radiation the corresponding weight is 

W z = l- z(l- z)[2- z(l- z)], (28) 

This weight is very well-behaved: ^ < W z < 1. The virtual-gluon contribution has to be 
compensated with the weight 



W A = exp(A G ), A G = $ 



G 



G 



(29) 
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3.3 Quark-gluon transitions 



Although the original CCFM equation was formulated for gluons only, we may try to 
extend it to quarks and allow for quark-gluon transitions. The treatment of the quark 
kernels is identical to the case of gluon kernels described in the previous subsections. As 
for the case of QG transitions, we intend to apply the same importance sampling for 
generation of the t and z variables as described in Subsection 13.11 using at the primary 
MC generation level the following approximate kernels 



ai 0) (t + ln(l-z)) e(l-z-e(t)) 



z?? K (z, t; t ) = K - ' ^ ~" ^ {_ z ^" [SikAP k + F 7 ( 2j , (30) 

where F^ = ma.x z Fj^ (z) are given in Appendix C of Ref . [23] . This means that at the 
low MC level we artificially include singular factors 1/(1 — z) for non-diagonal transitions 
that are not present in the corresponding exact kernels. The above approximation is then 
compensated by the MC weight being the ratio of the exact to approximate kernels 

The loss of efficiency due to this artificial modification is rather small, whereas the gain 
in simplicity of the algorithm is significant. 

For the exponent of the Sudakov form-factor we get 

2 i ,o, 



*i! *o) = ^ { a kk [9(h) - g(h) 



(32) 




where the functions Fj^-(z) are given Appendix C of Ref. [23]. Again, the non-singular 
terms in the above form factor have to be integrated numerically. 



3.4 Numerical tests 

We have performed comparisons of the MC solution of the above evolution equations, im- 
plemented the program EvolFMC [32], with the solution provided by the non-MC program 
APCheb40 [33] in which the same ^-dependent strong coupling has also been implemented. 
In both cases we have evolved the singlet PDFs for gluons and three doublets of mass- 
less quarks from Q = 1 GeV to Q = 10, 100, 1000 GeV. We have used the following 
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parameterization of the starting parton distributions in the proton at Qo — 1 GeV: 



xD G 




= 1.9083594473 


. x -°- 2 (l-xf-° 


xD q 




= 0.5 ■ xD sea (x) 


+ xD 2u (x), 


XDg 




= 0.5 • xD sca (x) 


+ xD d (x), 






= 0.6733449216 


. x -°- 2 (l-x) 7 - 


xD 2u 




= 2.1875000000 


■x°- 5 (l-x) 3 - , 


xD d 


\x) 


= 1.2304687500 


■x°- 5 (l-x) 4 - , 



(33) 





Figure 1: The upper plot shows the singlet quark distribution xD q (x,Q) evolved from Qo = 
lGeV (black) to Q = 10 (red), 100 (green) and 1000 (blue) GeV, obtained from EvolFMC 
(solid lines) and APCheb40 (dashed lines, hardly distinguishable), while the lower plot shows 
their ratio. 

In Fig. [TJ we show the resulting quark distributions obtained from the two programs 
(the upper plot) as well as their ratios (the lower plot). As one can see, these two 
calculations agree at the level of 0.2%, except for the x-region very close to 1, where 
the MC statistics is low. The similar agreement has been found for the resulting gluon 



9 




Figure 2: The upper plot shows the gluon distribution xD q (x,Q) evolved from Q = 1 GeV 
(black) to Q = 10 (red), 100 (green) and 1000 (blue) GeV, obtained from EvolFMC (solid 
lines) and APCheb40 (dashed lines, hardly distinguishable), while the lower plot shows their 
ratio. 

distributions, shown in Fig. [2l The precision of the presented results is limited by the 
statistical errors of the MC calculations. 



4 One-loop CCFM equations of Kwiecinski 

The CCFM equations [12] describe the evolution of an unintegrated parton distributions 
cIk(x, k T ,Q) which depend on the parton transverse momentum k T in addition to the 
longitudinal momentum fraction x and the scale Q. They are related to the integrated 
parton distribution of the DGLAP equations Dk{x,Q) through the following relation 

xD K (x,Q) = J d 2 k T xd(x,k T ,Q) (34) 

Originally, the CCFM equations were derived for the unintegrated gluon distribution. 
In the all-loop approximation, they also take into account angular ordering (coherence) 
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in the initial state gluon cascade for both large and small values of x. This leads to a new 
non-Sudakov form factor that sums virtual corrections for small x. 

In the one-loop approximation, the angular ordering at small x and the corresponding 
virtual corrections are subleading. Thus, the resulting equations contain only coherence 
in the real parton emissions for large x, including only the Sudakov form factor which 
resumes virtual corrections. At small x the standard DGLAP transverse momentum 
ordering appears. The one-loop CCFM equations for both quark and gluon unintegrated 
distributions, written for the first time in [15], take the following form 



d 2 q 



9(q - qo)0(Q - q) J y z? KJ (q, z) djf^, l T , q 

X 

P = k T + (1 - z )q. (35) 

This equation is slightly more general than the original equation from Ref. [15] since 
it allows for the coupling constant (hidden in the splitting functions T^j(g, z)) which 
depends both on z and q. The iterative solution to eq. fl35|) reads 

i 

xdx(x, k T , Q) = xd° K (x, k T ) + / dz / d 2 k^ 

n=l { J Ko,...,^-! 



II J 2^^ qi ~ Qi ~^ J dz i Z ^ K.K^Mii Z i) xd° Ko (^Z ,kQ 
n n 

Sfx - J]*) 5 2 {k T + ^(1 - Zi )q>) , (36) 



i=0 i=l 



where K = K n . As compared to the basic iterative solution of the DGLAP-type equation, 
eq. (T5]), the above series differs only by the presence of the independent angular integrals 
in d 2 qi and by the presence of the delta function 5 2 {k T — . . . ), which in the Markovian-type 
algorithm plays only a "spectator" role. In addition, in eq. f)36p the Sudakov form factor 
has not been explicitly worked out, as in eq. ([2]), and the initial condition xd° K is now 
/c T -dependent. Having that in mind, we reorganize eq. (136]) in the form similar to the 



11 



DGLAP solution, eq. (J2J, 

1 



xd K (x, k T , t) = e- 9K ^x(f K (x, k T , t ) + ^ / rfxo / ci 2 ^ ^ 

n=l ^ ^ K ,...,K n -x 

n * 1 n 

II / ^6(^-^0 J^Jdz, e-^^H Zi)e-*^«^ 



to 



X 



n n 

xd K (^o, , ^o) S(x - x Y\ Zi) 5 2 (k T - kl + ^(1 - Zi)q t 



i=l i=l 



(37) 



where as before % = exp(tj) and 0j is the azimuthal angle of Notice that t = InQ 
plays now the role of the evolution variable. 

It is now transparent that as described in Ref. [23] the forward Markovian algorithm 
EvolFMC can easily be extended to embed the generation of parton transverse momentum 
k T . After n steps of the forward Markovian evolution the transverse momentum of the 
off-shell parton entering the hard process becomes 

n 
i=l 

where k$ is an intrinsic parton transverse momentum and cfi a 2-dimensional evolution 
variable. The physical transverse momenta of emitted particles are 

(39) 

In the MC evolution q is constructed as follows: 

q { = e tl (cos 0;, sin 0;), (40) 

where ti is an evolution variable and fa is an azimuthal angle generated at each evolution 
step from a flat distribution in the range [0, 2ir}. 

In our numerical tests, which will be reported in the next section, the intrinsic parton 
transverse momentum k$ is generated at the initial evolution scale to from the following 
x-independent distribution 



i-^h (41) 

where k is some adjustable parameter, set to 1 GeV in our tests. 



5 Numerical results for /^-distributions 

In the following we present results for the unintegrated parton distributions as func- 
tions of the transverse momentum k T , generated by the forward Markovian MC program 
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EvolFMC [32] according to the one-loop CCFM equation described in the previous section. 
The starting point of the evolution is Qo = 1 GeV and the initial conditions are specified 
in eqs. (1531 and (jUp. 

In Figs. [3] and H] we show the gluon distribution obtained from the CCFM equation 
with gluons only, as given by Marchesini and Webber in Ref. [13]. The distribution was 
additionally integrated over x. For rising values of the scale Q such a gluon distribution 
moves towards large values of k T , becoming at the same time less sensitive to the low 
k T region. This observation is summarized in Fig. [5] where we present the average gluon 
(k T ) 2 as a function of x for Q = 1, 10, 100 and 1000 GeV. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 

log 10 (k T [GeV]) 



Figure 3: Gluon k T distributions integrated over x for the one-loop CCFM equation of 
Marchesini-Webber, obtained from EvolFMC for Q = 1 (black), 10 (red), 100 (green) and 
1000 (blue) GeV. 

In Figs. M and El we show the results of the same studies for the full one-loop CCFM 
equation for the singlet quark and gluon distributions, given by Kwiecihski et al. in 
Ref. [15]. The change of the distributions (integrated over x) with the hard scale Q is the 
same as in the previous case. Both the quark and gluon distributions become less sensitive 
to the "soft" k T values moving towards the region of large k T, s. This is also visible in 
Fig. M for the average transverse momentum for gluons and quarks. As expected, the 
average k T rises when x — > [15]. 

6 Summary and outlook 

In this work we have presented a systematic extension of the Monte Carlo algorithm that 
solves the DGLAP equation into the algorithm solving the one-loop CCFM equation. 
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Figure 4: Gluon k T distributions multiplied by (k T ) 2 integrated over x for the one-loop CCFM 
equation of Marchesini-Webber, obtained from EvolFMC for Q = 1 (black), 10 (red), 100 
(green) and 1000 (blue) GeV. 



Gluons <kj> 




-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

iog 10 (x) 



Figure 5: The average (k T ) 2 of gluon as a function of x for the one-loop CCFM equation 
of Marchesini-Webber, obtained from EvolFMC for Q = 1 (black), 10 (red), 100 (green) and 
1000 (blue) GeV. 



To this end two technical problems have been solved: the coupling constant has become 



14 



dx xD G (x,k T ,Q) 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 
"og 10 (k T [GeV]) 
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dx xD a (x,k T ,Q) 




' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 1 1 ' ' ' I 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 
log 10 (k T [GeV]) 



Figure 6: Gluons and quarks k T distributions integrated over x for the one-loop CCFM 
equation of Kwiecihski, obtained from EvolFMC for Q = 1 (black), 10 (red), 100 (green) and 
1000 (blue) GeV. 




k T [GeV] k T [GeV] 

Figure 7: Gluons and quarks k T distributions multiplied by (k T ) 2 integrated over x for the 
one-loop CCFM equation of Kwiecihski, obtained from EvolFMC for Q = 1 (black), 10 (red), 
100 (green) and 1000 (blue) GeV. 



^-dependent and the evolution in transverse momenta has been added, in addition to 
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Gluons <ky> Quarks <kj> 




log 10 (x) log 10 (x) 



Figure 8: The average (k T ) 2 of gluons and quarks as a function of x for the one-loop CCFM 
equation of Kwiecihski, obtained from EvolFMC for Q = 1 (black), 10 (red), 100 (green) and 
1000 (blue) GeV. 

the evolution in the longitudinal momentum fractions x in the parton distributions. The 
modification of the coupling constant has lead to the need of solving numerically the 
transcendental equation. First numerical results have been presented, confirming the 
known observations on CCFM unintegrated PDFs. 

The presented algorithm for the one-loop CCFM evolution we consider as the first 
step in extending the other types of the QCD evolution equations available in our MC 
programs beyond that of the DGLAP type. In particular, implementing the complete 
CCFM evolution, both in the Markovian and non-Markovian (constrained) algorithms, 
see ref. [26], is now in an advanced stage. The presently implemented one-loop CCFM 
option will be used in the forthcoming studies of various aspects of the QCD evolution 
equations according to several different schemes. 
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